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LECTURE 3. FIRST-ORDER LINEAR EQUATIONS 



First-order linear differential equations. We will give a systematic method of solving first-order 
differential equations (of normal form) 

(3.1) y' +p{x)y = f{x) 

on a given interval I, where p, f are continuous functions. 
First, the homogeneous equation 

(3.2) y'+p{x)y = 

is solved by quadrature. Let P{x) = f p{x)dx be an indefinite integral oip{x). Then, 

-l(e^Wy) = e^Wy' + p{x)e^^^^y = e^^^\y' + p{x)y) = 
ax 

if and only if y solves (3.2), since e^^^^ / 0. 

Theorem 3.1. Let p{x) be continuous on an interval I and P{x) = / p{x)dx be an anti-derivative ofp{x). 
Then, 

is a solution of {3.2) for any constant c. Conversely, all solutions of (3.2) are of this form. 

Exercise. Show that a solution y{x) of (3.2), where p{x) is continuous, is: either y{x) = for all x 
or y{x) / for all x. 

Next, we treat the inhomogeneous equation 

(3.3) y' +p{x)y = f{x) 

by variation of parameters. Let P[x) = J p{x)dx as before. Then, 

^(e^(^)y) = e^(^)(y' + p(x)y) = e^(^)/(x). 

/•X 

Hence, ye^^^^ = Vo + e^^^^f{s)ds for some yo and for some xq. 

JXQ 

Theorem 3.2. Let p{x) be as in Theorem 3.1 and let f{x) be continuous on an interval xq g /. Then, the 
general solution of (3.3) zs given by 

y{x) = yoe"^(^) + e-^(^) T e^^''^ f{s)ds. 

Jxo 

Moreover, y{xo) = yo if and only if P{x) = /^^^ p{s)ds. 
Example 3.3. Consider the differential equation 

(3.4) y' +y = x + 2,. 

By trying y{x) = ax + h one easily finds a solution y = x + 2 of (3.4). If y = (t){x) is another solution 
of (3.4), then z = y — {x + 2) must satisfy the corresponding homogeneous equation z' + z = 0. By 
Theorem 3.2, then, z = ce~^. Therefor, the general solution of (3.4) is 

y = ce~'-^ + X + 2. 
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Exercise. (Bernoulli equations) 1. If p and q are continuous functions of x and n 7^ 1 is constant, 
show that the Bernoulli equation 

y' + p{x)y = q{x)y" , y>0 

can be reduced to a linear equation upon substitution y = u"^ for a suitable constant m. The 
condition y > ensures that the function u = y^/™ is meaningful. When n = 1 the original 
equation is already linear. 

2. Solve the differential equation y' + y = xy^ , y > 0. 



The logarithmic spiral. Suppose a curve r = f{d) in polar coordinates cuts the radius at a con- 
stant angle, say ip, as shown in the figure below. 




Figure 3.1. The logarithmic spiral. 



The game is to find the equation of this curve. If ip = the curve is a ray extending from the 
origin to 00 and it cannot be represented in the form r = f{9). If tp = 7r/2 or — 7r/2 the curve is a 
circle centered at he origin. If 

(3.5) -7r/2 < V < vr/2, iP^O, 



we write tan ip = —. The differential equation is 



tan ip 



(See Figure 3.2.) 
Using tan = 

(3.6) 



we obatin 



ke 



rdO 
dr 



1 



rd9 
dr 



dr 
d9 




kr. 



Figure 3.2. 



and therefore r = ce 

A curve of this kind is called a logarithmic spiral. The steps in the discussion above are reversible, 
and hence a curve r = f{0) cuts the radius at a constant angle -0 satisfying (3.5) if, and only if, it 
is a logarithmic spiral. Since k is arbitrary, this gives a geometric interpretation of all processes 
obeying the exponential law of growth. 

A logarithmic spiral looks like a snail shell (Figure 3.1), and this is not a coincidence. A snail 
shell has the characteristic that it grows only at one end, and the part of the shell already laid 
down does not change. The growth is specified by its rate in two perpendicular directions - radial 
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and transverse. Let us suppose that both rates are proportional to the size at time t. The precise 
meaning of "size" will not matter, and we use the weight as a quantity that is easily measured. 
If arc length in the radial direction is si and in the transverse direction is S2, then 

(3.7) ^ = kiW and ^ = W, 

at at 

where ki , ^2 are constants. The arc length in polar coordinates thus satisfies 

ds^ = dr^ + {rdef. 

Setting dO = we get dsi = dr and setting = we get ds2 = rdO. Dividing the second equation 
in (3.7) by the first, we arrive at 

rdO /c2 

dr ki 

This is the same equation as in (3.6) with k = ki/k2. 
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